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^ ■ Abstract 

' We extend previous relativistic mean-field (RMF) calculations of multi-X nuclei, using vector 

boson fields with SU(3) PPV coupling constants and scalar boson fields constrained phenomeno- 
logically. For a given core nucleus, the resulting K separation energy Bj^, as well as the associated 
nuclear and X-meson densities, saturate with the number n oi K mesons for k > Kgat ~ 10. Sat- 
uration appears robust against a wide range of variations, including the RMF nuclear model used 
^ , and the type of boson fields mediating the strong interactions. Because Bj^ generally does not 

pH ' exceed 200 MeV, it is argued that multi--fr nuclei do not compete with multihyperonic nuclei in 
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providing the ground state of strange hadronic configurations, and that kaon condensation is un- 
likely to occur in strong-interaction self-bound strange hadronic matter. Last, we explore possibly 
self-bound strange systems made of neutrons and mesons, or protons and K~ mesons, and 
study their properties. 
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I. INTRODUCTION AND OVERVIEW 



Kaon condensation in dense matter was proposed over 20 years ago by Kaplan and Nelson 
. It is necessary to distinguish in this context between K mesons and K mesons which 
interact quite differently with matter. The empirical evidence from K~ atoms is that the 
A'-nuclear interaction is strongly attractive, and absorptive as well, with typical values of 
150 — 200 MeV attraction at nuclear-matter density po, as reviewed recently by Friedman and 
Gal [3]. A strong nuclear attraction of somewhat less than 100 MeV at po for K~ mesons, 
compared to a weak repulsion of order 25 MeV for mesons, follows from observations of 
enhanced near-threshold production of K~ mesons in proton- nucleus collisions at GSI {4]. 
This weakly repulsive nature of the i^'^-nuclear interactions was quantified a long time ago, 
starting with Dover and Moffa j^, and is also reviewed in Ref. ^. Given the distinction 
between the nuclear interactions of K mesons and K mesons, the term kaon condensation 
is used loosely here and elsewhere to mean K condensation. 

Neutron stars, with a density range extending over several times ppjoffer the most natural 

dense systems where kaon condensation could be realized; see Refs. HflflSQ,! [ll| for com- 
prehensive reviews of past work. We note that in Heaven, for neutron stars, weak-interaction 
time scales of order 10"^ s and longer are operative, enabling strangeness-changing processes 
such as + v^. to transform high-pressure dense electrons to K~ mesons once the 

effective mass of K~ mesons drops down below approximately 200 MeV. Under some opti- 
mal conditions, recalling that K mesons undergo attraction of order 100 MeV per density 
unit of Po kaon condensation could occur at densities about 3po, depending on the way 
hyperons enter the constituency of neutron stars as first recognised by Ellis, Knorren and 
Prakash [isl. However, on Earth under laboratory conditions, strong-interaction time scales 
of order 10~^^ s are operative; processes of equilibration and hadronization subsequent to 
dense-matter formation in heavy-ion collisions occur over much shorter times than those 
controlling the composition of neutron stars. If antikaons bind strongly to nuclei, according 
to a scenario spelled out recently by Yamazaki et al. [l4|, then K mesons might provide 
the relevant physical degrees of freedom for self-bound strange hadronic matter that would 
then be realized as multi-^ nuclei. It requires that the K separation energy B^^ beyond 
some threshold value of strangeness exceeds uikc^ + — m/^c^ > 320 MeV, where Pat is 
the nucleon chemical potential, thus allowing for the conversion A ^ ^ + in matter. 
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For this strong binding, A and S hyperons would no longer combine macroscopically with 
nucleons to compose the more conventional kaon-free form of strange hadronic matter fisl . 
K mesons will then condense macroscopically. However, the nuclear densities encountered 
in these strange hadronic nuclei are somewhat less than the typical 3po threshold required to 
lower sufficiently the K energy in matter to reach condensation. Yet, precursor phenomena 
to kaon condensation in nuclear matter could occur at lower densities as soon as Bf^ exceeds 
the combination mK(? + ^in — ^t.(? ~ 240 MeV. In this case, the only mechanism underlying 
the widths of multi-^ states is the fairly weak conversion KNN AN. 

Recently we have reported on preliminary calculations of multi-^ nuclear configurations 
[igI using the relativistic mean-field (RMF) methodology, constrained by ^-nucleus phe- 
nomenology. It was found that the nuclear and K densities behave regularly on increasing 
the number of antikaons embedded in the nuclear medium, without any indication for abrupt 
or substantial increase, and that the K separation energy saturates. Roughly speaking, the 
heavier the nucleus is, the more antikaons it takes to saturate the separation energies, but 
even for ^osp]^ ^^le number required does not exceed approximately 10. Because the calcu- 
lated K separation energies Bf^ do not generally exceed 200 MeV, for input binding in the 
accepted "deep-binding" range -B^- ~ 100 — 150 MeV for a single K meson [l7, 1^, 19|, it was 



deemed unlikely that kaon condensation occurs in nuclear matter. This leaves antikaons in 



multi-i^' nuclei com: 
hadronic matter 



brtably above the range of energies appropriate to (hyperonic) strange 
15[ |. In the present article we discuss the full scope of these calculations 
demonstrating the robustness of this saturation property. In particular we study the sen- 
sitivity of the results to the nuclear equation of state used, through the nonlinear RMF 
version employed, and the role of "hidden strangeness" isoscalar meson fields beyond the 
standard isoscalar, scalar (a), and vector [uj) meson fields. Although both cr- and cu-meson 
fields mediate attraction between K mesons and nucleons, they play different roles for the 
interactions within K mesons, similarly to the pattern well known for nucleons. The a meson 
induces attraction, whereas the u meson induces repulsion. If the ^-meson couplings were 
exclusively limited to scalar-meson fields, the resulting ^-meson separation energies would 
not have saturated. However, chiral model studies of KN low-energy phenomenology give a 
clear evidence in favor of the lowest-order Tomozawa- Weinberg vector interaction, which in 
terms of meson exchanges is equivalent to vector-meson exchanges with purely F-type SU(3) 
pseudoscalar-pseudoscalar- vector (PPV) vertices [3]. Our philosophy in this work is to use 
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these vector-meson fields coupling constants as they are, augmenting the ^-nucleus vector 
interaction by additional scalar couplings such that Bj^ ~ 100 — 150 MeV holds for single-^ 
nuclei. We find no precursor behavior to kaon condensation for K mesons in self-bound 
nuclear matter. 

We also explore in this work exotic strange self-bound configurations where K mesons are 
bound to either neutrons or protons. The simple example of a quasibound K~pp system (and 
thus also its charg e-sy mmetric partner K^nn) recently calculated solving Faddeev coupled- 



channel equations |20l . l2ll . |22| | , clearly demonstrates that K mesons can bind together nuclear 
clusters that are otherwise unbound. The point here is that the underlying K~p and K^n 
interactions (each with equally mixed / = and 1=1 components) provide considerably 
more attraction than the purely 1=1 K^n and K^p interactions. The RMF calculations 
reported here start with eight neutrons, showing that a finite number of neutrons can be 
made self bound by adding together a few mesons, with K separation energies of order 
Bj^ ~ 50 — 100 MeV. We study the role of the isovector p meson in stabilizing these exotic 
configurations, owing to its role in distinguishing between the underlying / = and 1=1 
KN interactions. We find that the emergent stable neutron configurations are more tightly 
bound than in the corresponding ordinary nuclei with N Z along the stability valley, and 
the neutron single-particle spectra display substantial rearrangement. However, these exotic 
configurations are found to be unstable against charge-exchange + n — > K~ + p 
reactions. 

In Sec. [IT] we briefly outline the RMF methodology and discuss the K coupling constants 
to the meson fields used in the present calculations. Results are shown and discussed in 
Sec. Unifor K separation energies and density distributions, also displaying the dependence 
on the type of nonlinear RMF model used and the contribution of specific meson fields to the 
energy systematics and particularly to maintaining saturation in a robust way. A separate 
subsection is devoted to the study of exotic multi--S' "nuclei" with neutrons only. Again, 
binding energies and densities are discussed, plus rearrangement features of the neutrons in 
the ^-extended mean field. We conclude with a brief summary in Sec. IIVI 
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II. METHODOLOGY 



A. RMF equations of motion 

Bound nuclear systems of nucleons and several K mesons are treated in this work within 
the RMF framework, where the interactions among hadrons are mediated by the exchange 
of scalar- and vector-meson fields. The model Lagrangian consists of a standard nuclear 
part Cn and the Lagrangian density Ck describing the kaonic sector: 

Ck = (V^Xy (V ^^K) - mlK^K + g^KfiiKK^K a + g^.K^iKK^K a\ (1) 

Here, 

1)^ = 8^ + 1 g^K uj^, + igpKr ■ pf, + ig^K(t)i, + ''^e\{l + rs)^^ , (2) 

and K {K'^) denotes the kaon (antikaon) doublet. To be specific, we discuss ii'~-nuclear 
systems. Similar expressions hold for mesons. In addition to a scalar-meson field a and 
to vector-meson fields uj and p normally used in purely nuclear RMF calculations, we also 
considered meson fields that couple exclusively to strangeness degrees of freedom, a scalar 
cr*, and a vector cj). Standard techniques yield a coupled system of equations of motion for 
nucleons and all meson mean fields involved; we refer the reader to our earlier work [16] for 
details. Here it suffices to recall that the presence of K meson(s) induces additional source 
terms in the Klein-Gordon (KG) equations for the meson (mean) fields. In the case of K~ 
mesons, the source terms contain the K~ density 

pK- = ^{Ek- + g^K t^o + QpK Po + g<f>K + e Ao)K~K'^ , j d?x pK- = , (3) 

where E^- = idtK~ . Hence, the K mesons modify the scalar and vector potentials that 
enter the Dirac equation for nucleons, thus leading to rearrangement of the nuclear core. The 
polarized nucleons, in turn, modify the ^-nucleus interaction. This calls for a self-consistent 
procedure for solving the equations of motion. 

In our model, the KG equation of motion for the K~ meson acquires the form 

[- - El- + + Re YIk- = , (4) 
with the K~ self-energy given by 

Re Hk- = - {gaKrriKO-o + ga*KmK(Tl) - 2Ek- {gu^K^o + 9pKPo + fl'</>A'0o + eAo) (5) 
- {guK^^o + gpKPo + 9^k4>o + eAof . 
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Of the three terms on the right-hand side (rhs) of Eq. ([5]), the first one is a scalar-meson 
contribution, whereas the other two terms are vector-meson contributions. The scalar con- 
tribution is sometimes lumped together with the kaon mass niK to form a density-dependent 
effective kaon mass via 

m*f^'^ = niK^ - gaKrriKcro - ga^KniKC^l . (6) 

Finally, to account for K~ absorption in the nuclear medium, the self-energy 11^^- = 
2EK'V-^yi-p i'^ Eq. (jl]) was made complex by adding IthYIk- and the real energy Ek- was 
replaced by E^- — iV k- /'i- The imaginary part of the self-energy, ImllR--, was taken from 



optical model phenomenology, with a strength fitted to atomic data j23j and with en- 
ergy dependence that follows the reduced phase space for the decaying initial state. We 
assumed two-body final-state kinematics for the decay products in the absorption channels 
KN vrF {Y = E, A) (80%) and KNN YN (20%) with branching ratios indicated in 



parentheses 
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The set of coupled equations containing the Dirac equation for nucleons, the KG equa- 
tions for the meson mean fields and for antikaons was solved fully self-consistently using an 
iterative procedure. This appeared crucial for the proper evaluation of the dynamical effects 
in nuclei with k {k, = 1,2,3, ■■■ ) K mesons. The K separation energy 

Bk = B[A,Z,kK]- B[A,Z,{k-1)K] , (7) 

where B{A, Z, nK) is the binding energy of the K^-nuclear system, contains mean-field 
contributions due to rearrangement of the nuclear core. 



B. Choice of parameters 



For the nucleonic Lagrangian density /Zat we used the RMF parameter sets NL-SH 24 1 



and NL-TM1(2) [25| which have been successfully used in numerous calculations of various 
nuclear systems. For the (anti-)kaon coupling constants to the vector-meson fields, we used 
a purely F-type SU(3) symmetry, ay = E/{E + D) = 1: 

2gu;K = V2g^K = 2gpK= gpTv = 6.04 , (8) 

where the value of (^p^ is due to the p ^ 2tt decay width. As mentioned in Sec. [U this choice 
corresponds to the underlying Tomozawa- Weinberg lowest-order term in chiral perturbation 



theory 19|. The value of g^^^ = 3.02 adopted here is considerably lower than the quark- 
model (QM) value applied to NL-SH, g^x = \g^N~^^ = 4.32, which was used in our previous 



work 



16| . and we consider it to be the minimal value suggested by theory. The K RMF 



vector potential at threshold in nuclear matter is then given, in the static approximation, 
by using the last two terms on the rhs of Eq. ([5]) : 

SU(3) NL-SH ^ . /nSU(3)„NL-SH^ \2 

T/K- _ iiuK {JuiN HO -L / y^j^ y^pj po \ _ ^ i^A^\r fn\ 

ml 2mK \ ml J 

with m^ = 783 MeV and Po^"^" = 0.146 fm'^ We point out that the value g^^'^^ = 12-95 
is not far away from the value g^fp^'^ = 11.06 from the latest A^A^-potential fit by the 
Nijmegen group |26]. This latter value was obtained in that A^A^ analysis after allowing for 
part of the isoscalar vector-meson field strength to result from a combined p-vr exchange. We 
also studied the role of isovector K nucleus interactions by comparing the results of using the 
present SU(3) choice g^]^^^^ = 3.02 with results applying a QM universal isospin coupling 
to NL-SH: gpK = g^N~^^ ~ 4.38. This value is substantially higher than the Nijmegen 
potential fit value g^^^^"^ — 2.77, apparently to compensate for disregarding the almost four 
times higher value of the tensor coupling constant /^*^°^. 

SU(3) symmetry is not much of help in fixing the (anti-)kaon coupling to the scalar- 
meson field a, simply because the microscopic origin of the a field and its various couplings 
are not unambiguous. It has been shown recently that interpreting the a field in terms of 
a {J'", I) = (0"^,0) resonance in the utt-KK coupled-channel system leads to a vanishing 
KN forward-scattering amplitude at threshold, thus suggesting a vanishing contribution to 
the corresponding ^-nucleus optical potential j27|. However, even for the empirically large 
value of gcrN obtained in the A^A^ case ((7^^^°^ = 10.17) and also within the RMF description 
of nuclei, there is no consensus on its microscopic origin, except that QCD sum-rules do 
produce strong scalar condensates. Modern A^A^ potentials using chiral perturbation theory 
guidelines obtain a rather strong isoscalar-scalar two-pion exchange contribution involving 
excitation of A(1232) in intermediate states 28|]. A similar two-pion exchange contribution 
for ^A^, involving the excitation of K*{892), cannot be excluded at present. In the absence 
of QCD sum-rule determinations of gax, one relies for an order of magnitude estimate 
on simplified models such as the QM, giving rise to g'^^ = \gaN, which for the NL-SH 
model gives g^^^ = 3.48. The associated RMF K~ nuclear scalar potential, in the static 
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approximation, is given by: 

QM NL-SH s 

^R^^F-scaia. = "^^^^^ = -66-3 MeV , (10) 

using the values mg. = 526.1 MeV from NL-SH and Pq ~ O.Qpo, where Pq is the scalar 
density. Our choice of g^K is conceptually different, fitting g^^x to several selected values of 
K separation energy Bj^ in nuclear systems with a single K meson. These fitted values, all 
of which were considerably lower than the QM value, are specified in the next section. Thus, 
our scalar potentials are viewed as a supplement to the minimal vector potentials discussed 
above to scan over K nuclear binding energies in a given energy range, without imparting 
any microscopic meaning to these scalar potentials. We also tested the effect of adding 
another scalar-meson field that couples exclusively to strangeness, "hidden strangeness" a* 
meson with mass m^* = 980 MeV, and coupling constant ga*K = 2.65 determined from 



/o(980) ^ K+K- decay (l2| 



III. RESULTS AND DISCUSSION 



A. Saturation of K binding energies and hadronic densities 



Following the observation made in Ref. 161] that K binding energies, as well as nuclear 



and K densities, saturate on increasing the number n of K mesons, we have explored how 
robust this saturation is. In particular we studied, for several selected nuclei across the 
periodic table, the role of various components of the ^-nucleus interaction in establishing 
saturation and the sensitivity to the choice of the RMF model. Representative examples are 
shown in Figs. [H and O 

Figure [U presents the Is K~ separation energy B^- in multi-fC" nuclei ^^O + kK' as a 
function of the number k, of K~ mesons, using the NL-SH RMF parametrization, for several 
mean-field compositions of the K~ self-energy Eq. ([5]) with K~ vector-meson couplings 
given by Eq. ([8]). The upper group of curves is based on a value of g„K = 2.433 ensuring 
Bx- = 100 MeV for k, = 1. The (j),p,a* meson fields do not practically contribute in this 
case, whereas the Coulomb field adds a few MeV attraction and Iml^pt adds a few MeV 
repulsion. For k > 1, the various curves of the upper group diverge from each other: with 
respect to a "minimal" a+u model (open circles), the main contributors are the repulsive 
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FIG. 1: Is separation energy B^- in ^^0+kK~ as a function of the number k of antikaons 
in several calculations detailed in the text, as listed in the inset, using the NL-SH RMF nuclear 
model with go-K = 2.433 for the upper group of curves and gaK = 1.703 for the lower group. 



and p vector mesons, as judged by the curves marked by 
respectively. Given their contributions, which get larger with 
field, the a* meson field and ImVopt makes a small difference. 
ImVJ^pt makes a big difference for the lower group of curves 
both with g„K = 1.703 fitted to Bk~^ ^ 40 - 50 MeV for k = 
ImKpt magnifies its effect for relatively low values of B^- 
adding significant repulsion the lesser the value of B^- is. 
curve in Fig. [1]) leads to a rapid fall-off of B^- , terminating 



solid circles and open squares, 
K, the inclusion of the Coulomb 
However, the K~ absorptivity 
consisting of only two choices, 
= 1. The energy dependence of 
in the region < 100 MeV, 
This added repulsion (lowest 
the binding at k = 3, because 
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FIG. 2: Is separation energy Bj^- in ^^Csl+kK' , as a function of the number k of K~ 
mesons, calculated in the NL-SH (circles, solid lines), NL-TMl (triangles, dashed lines), and NL- 
TM2 (squares, dot-dashed lines) RMF models. The lower (upper) group of curves was constrained 
to produce B^- = 100 (130) MeV for k = 1. 

the system ^^0+4K~ is found to be unbound for this particular choice of k = 1 parameters. 
The lesson from Fig. [1] is that saturation of the binding energy in nuclear systems with 
K, i^'" mesons is a robust phenomenon, which remains valid regardless of the type of meson 
fields mediating the strong interaction among antikaons and nucleons, provided a minimal 
isoscalar vector-meson field (u) is included. For a sufficiently large number k of K" mesons, 
the combined repulsive K^K^ interaction generated by the vector meson fields u, 0, p wins 
over the attractive interaction generated by the isoscalar scalar-meson fields (dominated by 
a). The effect of adding the a* scalar field is found to be insignificant. These conclusions hold 
also for a Lagrangian in which scalar fields are introduced differently than in Eq. ([T]), resulting 
in a correspondingly different definition of effective masses, = uik — gaxo'o — ga*KO'o 



29j, than in Eq. 
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Figure [2] shows the Is K~ separation energy B^- in multi--ft'~ nuclei ^°Ca + nK~ as a 
function of the number k of K~ mesons, calculated in the NL-SH, NL-TMl, and NL-TM2 
RMF models for two choices of g^K designed, within each model, to produce B^- = 100 and 
130 MeV for k = 1. The values of g^K for NL-SH were 1.703 and 2.993, respectively. The 
difference between the various curves, for a given starting value of B^- , originates from the 
specific balance in each one of these RMF models between the vector fields and the scalar 
field. The figure illustrates that the saturation of the K binding energy in nuclear systems 
with several antikaons is not limited to a particular choice of RMF parametrization but is 
a general feature independent of the applied RMF model. Without loss of generality, we 
therefore specialize in the subsequent discussion to a specific RMF model, namely NL-SH. 

The dependence of the nuclear density pAr(r) and the K~ density PK~{f^) on the number 
K of K~ mesons in multi-i^'~ nuclei ^°Ca + k,K~ is shown in Fig. [31 The coupling constant 
gcrK = 1.703 was chosen such that the single- iiT" configuration was bound by 100 MeV, the 
same as for the NL-SH lower curve in Fig. [2l The density distribution pjv for "^°Ca is also 
shown, for comparison, by the dotted curve in the upper panel. It is clear from this figure 
that the central nuclear density pn saturates for k = 8 at a value about twice larger than 
that for in ^°Ca. In the lower panel, it is seen that the gradual increase of pk- with k 
slows down with increasing k. 

The saturation of the nuclear density in multi-^ nuclei manifests itself also in the behavior 
of the K effective mass in the nuclear medium, Eq. ([6]), as a function of the number k of 
antikaons. This is illustrated for ^°^Pb + k,K~ in Fig. HI Note that the calculated effective 
mass distribution m^_(r) remains almost independent of the number of K~ mesons over 
a large volume of the nucleus, for r > 3-4 fm, reflecting a similar k independence of the 
scalar a field through the underlying nuclear density. In fact, the a field in this region 
is almost unaffected by the presence of K~ mesons, as demonstrated by the dashed curve 
which uses the "static" ^°^Pb a field from a purely nuclear RMF calculation. It is only within 
a relatively small region near the nuclear center, typically r < 2-3 fm, that the variation 
of m^_(r) with n gets to be more pronounced. However, m*^^{r = 0) quickly saturates, 
already for k ^ 8. The figure demonstrates clearly that the concept of "nuclear matter" is 
far from being realized, even for a nucleus as large as ^°^Pb and that conclusions made on 
K binding and kaon condensation in finite nuclei, using nuclear-matter arguments, should 
be taken with a grain of salt. 
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FIG. 3: Nuclear density pN (top panel) and Is K density px- (bottom panel) in ^^Csl+kK , 
calculated in the NL-SH RMF model, with g^x = 1-703 chosen to yield B^- = 100 MeV in 
'^^Ca.+IK^ . The dotted curve stands for the '^''Ca density in the absence of mesons. 

B. Exotic K nuclear configurations 



Because in the underlying KN dynamics the / = interaction is considerably more 
attractive than the 1=1 interaction, we have looked for ways to maximize the role of 
the KN 1 = channel in multi-^ nuclei. For a nuclear core with N = Z, no matter 
which charge states are assigned to the K mesons, the average KN interaction is given 
by a (2/ + l)-average which disfavors the / = channel. This disadvantage is partly 
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FIG. 4: Is K~ effective mass m-^^ in mufti- nuclei '^^^Fh+KK~ , calcufated in the NL-SH model 
with g^K = 2.433 chosen to yield Bj^- = 100 MeV in ^'^^Pb+li^~. The dashed curve stands for 
the "static" case where the ^'^^Pb a fieid in a purefy nuclear RMF cafcufation was used in Eq. 
for m*^_. The dashed arrow indicates the charge half-density radius i?ch in ^'^'^Pb. 

removed by considering nK^ (or pK~) multi-^ nuclei, where both isospin channels assume 
equal weight, so that the stronger / = component may provide sufficient attraction to 
overcome the insufficient attraction in purely neutron matter. We therefore studied exotic 
configurations consisting solely of mesons bound to neutrons. Our calculations confirmed 
that K mesons can bind together systems of nucleons that otherwise are unbound. 

In Fig. [5l we compare the separation energies Bj^ in 16n+K,K^ and in 8n+K,K^ exotic 
multi-^ configurations with Bj^ in ^^0+k,K~ multi-^ nuclei, most of which were calculated 
in the NL-SH RMF model with the "canonical" Qvk coupling constants of Eq. ([8]) and 
g^K = 2.433 chosen to yield Bk- = 100 MeV in ^^O+IK' as in Fig. [H and for ImVopt = 0. 
For each sequence of multi-^ nuclei, Bj^ increases as a function of k to a maximum value 
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FIG. 5: Is K separation energy Bk in ^^0+kK-, 16n+KK°, and Su+kK'^, function of k, 

QT T/0\ 

calculated in the NL-SH RMF model, with gpK = (dot-dashed curves), g^j^ = 3.02 (solid 
curves) and QpK = QpN = 4.38 (dashed curves). ImVopt = is assumed everywhere except for the 
two lowest solid curves (open squares) in Sn+nK^ where ImVopt 0) such that the value of width 
F/fo is held fixed at 50 and 100 MeV, respectively, see text. 

and then starts to decrease. Whereas the sequence consisting of 8 neutrons plus mesons 
starts with k = 1 (not shown in the figure), a larger number of neutrons generally requires a 
threshold value for k, as shown for the sequences consisting of 16 neutrons plus K'^ mesons. 
Exceptions to the use of the canonical g^x set of Eq. ([8]), or to ImVopt = 0, are as follows: 

• 9^K^^^ — 3-02 was used everywhere except for Qpx = in the dot-dashed curves and 
except for gpx = QpN = 4.38 (universal isospin coupling) in the dashed curves to study 
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the role of the p meson in "nonexotic" multi-^ nuclei {}^0 + nK~) and in "exotic" ones 
{IQn+nK^). In ^'^0+k.K~ ^ the values of B^- for a given value of k > 1 decrease as gpK 
is increased, as expected from the repulsive K^K^ isovector interaction. In contrast, 
the larger the value of gpK is, the larger is the value of Bj^o expected in IGn+nK^, 
because it is the p isovector interaction that distinguishes the more attractive / = 
component of the K^n interaction from the less attractive 1 = 1 component. Indeed, 
this holds for /t < 8 in the figure. However, for k > 8, the contribution of the repulsive 
K^K^ isovector interaction becomes substantial for the values of QpK 7^ used here; 
the Bf^o dashed curve for the universal p coupling heads down, crossing the -B^o solid 
curve corresponding to SU(3) p coupling. All in all, substantial binding in lQn-\-nK^ 
multi-^ nuclei is reached for these values of Qpx 7^ 0. 

• The effect of ImVopt on B^^ is relatively unimportant for B^^ > 100 MeV, where the 
dominant KN — > vrS decay channel is closed. The inclusion of ImVopt is found then 
to induce repulsion of less than 5 MeV. However, in 8n+K,K^ multi-^ nuclei, where 
B^o < 80 MeV, the effect of ImVopt becomes significant. An estimate of this effect 
is given by comparing the Bj^o curve for F = (solid squares) with the Bj^o curves 
using ImVopt 7^ (open squares) such that the value of Tj^o is held fixed at 50 and 100 
MeV. As expected, the larger input widths induce a stronger repulsion that lowers the 
calculated Bj^o values. Yet considerably lower values of Bf^o are obtained once the K'^ 
widths are included self-consistently in these dynamical calculations. 

The nucleon-density distribution pAr(r) and the ^-density distribution p^^ (r) are shown 
in Fig. inifor ^^0+8K^ , 16r;,+8^'', and 8n+8K^. We note that p^v and pj^ are normalized 
to the number of nucleons and number of antikaons, respectively. The K couplings were 
chosen such that the 1K~ configuration in ^^O is bound by 100 MeV, as in Fig. [51 For 
comparison, we also present the density distribution p^ for ^^O without K mesons. Owing 
to the substantial K density p^ in the nuclear center, the central nuclear density PAr(O) in all 
three systems with 8 K mesons is about 2-3 times larger than the central nuclear density po 
in ^^O for k, = 0. The situation is particularly pronounced in 8n+8K^, with the same central 
density PAr(O) as in the systems with 16 nucleons + 8K. Furthermore, the 8n+8K^ system 
is compressed substantially in comparison with the other multi-^ systems, judging by the 
radial extension of p^ and p^ in both panels of Fig. [61 The relatively high value Pk{0) ~ 5po 
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FIG. 6: Nuclear density pN (top panel) and Is K density pp^ (bottom panel) in "^^O+SK , 
l%n+SK^ and Sn+SK^^ calculated in the NL-SH RMF model, with the "canonical" q^k cou- 
pling constants of Eq. ^ and with g^x = 2.433 to yield Bk- = 100 MeV in ^"^O+IK" as in 
Fig. [5j The dotted curve stands for the ^^O density in the absence of K mesons. 



for this system does not introduce complications due to possible overlap between antikaons, 
because the mean-square radius of K~ is less than half of the corresponding quantity for 
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Given the compressed nuclear densities plotted in Fig. El we show in Fig. [7] the calculated 
neutron single-particle energy levels in ^^O, in ^'^O+nK" and in Sn+zt^*^ multi-^ nuclei 
for K = 4, 8. The K couplings are the same as in Figs. and El again chosen to ensure 
Bx- = 100 MeV in ^^0+lK~ . The lsi/2 and 1^3/2 levels undergo increasingly attractive 
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FIG. 7: Neutron sing le-particle spectra in ^"^O (center) ^'^0+4(8)i^" (right) and 8n+4(8)K° (left), 
calculated in the NL-SH RMF model with g^x = 2.433 chosen to yield B^- = 100 MeV in 

shifts on varying k in these multi-^ systems. Particularly strong is the downward shift 
of the lsi/2 level, by about 70 MeV in ^^0+87^- and by about 130 MeV in 8n+8K^. In 
contrast, the lpi/2 neutron level is pushed up by about 10 MeV in the ^^0+k,K~ systems as 
a result of a gradually increasing spin-orbit splitting which reaches 43 MeV for k = 8 (recall 
that it is 7 MeV for k, = using NL-SH). The lpi/2 neutron level is weakly bound in the 
exotic 8n+K,K'^ systems for 1 < k < 6, getting more bound with k as shown in the figure for 
these systems. The Ip spin-orbit splitting becomes as large as 56 MeV in the exotic 8^ + 8^" 
system which exhibits the largest single-particle level splittings in this figure. Here the lpi/2 
neutron level, too, undergoes attraction. 

It is worth mentioning that exotic multi-K° configurations that contain no protons lie 
high in the continuum of "nonexotic" multi-^ nuclei that are based on nuclear cores with 
protons and neutrons. Figure [8] shows the calculated total binding energy B[A, Z, nK], 
Eq. ([7]), assuming for simplicity gpK = 0, for three sequences of multi--ft" nuclei. To illustrate 
the relationship between "exotic" and "nonexotic" configurations, we take the 16n+8K^ 
configuration, specifically in its lowest isospin 1 = 4 state, and replace successively K'^+n 
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FIG. 8: Total binding energy B[A, Z, kK] of ^^0+kK (circles), 16n+KK^ (solid triangles) and 
8n + kK^ (squares) multi--^ systems, as a function of k, calculated in the NL-SH RMF model with 
g^j^ = 2.433 chosen to yield B^- = 100 MeV in and with gpx = 0. For k = 8, total 

binding energy values for configurations that are intermediate between ^^0+8K^ and 16n+8^'^ 
are shown in open triangles along the dashed line. 

pairs by K~+p pairs until ^^0+8K~ is reached. This is demonstrated by the empty triangles 
along the vertical dotted line that connects the initial and final configurations. Both initial 
and final configurations have identical quantum numbers B = 8, Q = 0, / = 4, so they 
are commensurate. Therefore, although K mesons are capable of stabilizing purely neutron 
configurations, these exotic configurations do not compete energetically with multi-^ nuclei 
based on nuclear cores along the nuclear valley of stability. 
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IV. SUMMARY 



In the main part of this work, we studied several dynamical aspects of multi-^ nuclear 
states within RMF methodology. In particular, we discussed in detail the saturation pat- 
tern of K separation energies and nuclear densities on increasing the number of antikaons 
embedded in the nuclear medium. Saturation was demonstrated to be a robust feature of 
multi-^ nuclei. The saturated values of Bj-^, for "natural" values of meson-field coupling 
constants were found generally to be below 200 MeV, considerably short of the threshold 
value a; 320 MeV needed for the onset of kaon condensation under laboratory conditions. 
We conclude, consistently with our earlier conjecture [igI, that K mesons do not provide 
the physical "strangeness" degrees of freedom for self-bound strange dense matter. 

We first explored contributions of specific meson mean fields to the K separation energy 
Bj^. Saturation of Bf^ emerged for any boson-field composition that includes the domi- 
nant vector cij-meson field, using the "minimal" SU(3) value suggested by the leading-order 
Tomozawa- Weinberg term of the meson-baryon effective Lagrangian. Moreover, the contri- 
bution of each one of the vector (/)-meson and p-meson fields was found to be substantially 
repulsive for systems with a large number of antikaons, reducing the K separation energy 
as well as lowering the threshold value of number of antikaons required for saturation to 
occur. In contrast, the Coulomb interaction and the addition of a hidden-strangeness scalar 
cT*-meson field have little effect on the binding energy balance and on the pattern of satura- 
tion. 

We also verified that the saturation behavior of Bj^ is qualitatively independent of the 
RMF model applied to the nucleonic sector. The onset of saturation was found to depend 
on the atomic number. Generally, the heavier the nucleus is, the more antikaons it takes to 
saturate their separation energies. 

The saturation phenomenon found for the K separation energy is also reflected in the 
nucleon and antikaon density distributions, with the assertions made above remaining valid. 
The saturation of the nuclear density in multi--R' nuclei manifests itself in the behavior of the 
K effective mass distribution in the nuclear medium. We stress that in the case of antikaons 
the concept of nuclear matter is far from being realized even in such a heavy nucleus as 
^°^Pb. Specifically, the reduction of m*^(r) on adding K mesons is pronounced only within 
a small region around the nuclear center. 
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In the second part of this work, we studied exotic configurations consisting exclusively of 
neutrons and X° mesons. We demonstrated that a finite number of neutrons can be made 
self-bound by adding few X° mesons, with the resulting nuclear configurations more tightly 
bound than ordinary nuclei. Saturation of Bj^o was found for these exotic configurations 
too. Yet, these exotic configurations consisting exclusively of neutrons and mesons lie 
high in the continuum of the less exotic multi-^ configurations based on nuclear cores along 
the nuclear valley of stabihty. 

In conclusion, over a wide range of variations our calculations do not indicate any pre- 
cursor phenomena to kaon condensation in self-bound strange nuclear systems. 

Acknowledgments 

This work was supported in part by the GA AVCR grant IAA100480617 and by the Israel 
Science Foundation grant 757/05. 



[1] D.B. Kaplan, A.E. Nelson, Phys. Lett. B 175, 57 (1986); ibidB179, 409 (1986). 
[2] A.E. Nelson, D.B. Kaplan, Phys. Lett. B 192, 193 (1987). 

[3] E. Friedman, A. Gal, Phys. Rep. 452, 89 (2007); and earlier references cited therein. 
[4] W. Scheinast et al, Phys. Rev. Lett. 96, 072301 (2006); and earlier references to nucleus- 
nucleus experiments cited therein. 
[5] C.B. Dover, P.J. Moffa, Phys. Rev. C 16, 1087 (1977). 
[6] C.-H. Lee, Phys. Rep. 275, 255 (1996). 

[7] M. Prakash, I. Bombaci, M. Prakash, P.J. Ellis, J.M. Lattimer, R. Knorren, Phys. Rep. 280, 
1 (1997). 

[8] H. Heiselberg, M. Hjorth- Jensen, Phys. Rep. 328, 237 (2000). 
[9] H. Heiselberg, V.R. Pandharipandc, Annu. Rev. Nucl. Part. Sci. 50, 481 (2000). 
[10] N.K. Glendenning, Phys. Rep. 342, 393 (2001). 

[11] A. Ramos, J. Schaffner-Bielich, J. Wambach, Lect. Notes Phys. 578, 175 (2001). 

[12] J. Schaffner, I.N. Mishustin, Phys. Rev. C 53, 1416 (1996). 

[13] P.J. Ellis, R. Knorren, M. Prakash, Phys. Lett. B 349, 11 (1995). 



20 



[14] T. Yamazaki, A. Dote, Y. Akaishi, Phys. Lett. B 587, 167 (2004). 

[15] J. Schaffner-Bielich, A. Gal, Phys. Rev. C 62, 034311 (2000); and references cited therein. 
[16] D. Gazda, E. Friedman, A. Gal, J. Mares, Phys. Rev. C 76, 055204 (2007). 
[17] J. Mares, E. Friedman, A. Gal, Phys. Lett. B 606, 295 (2005). 
[18] J. Mares, E. Friedman, A. Gal, Nucl. Phys. A 770, 84 (2006). 

[19] W. Weise, Proc. IX Int. Conf. Hypernuclear and Strange Particle Physics, Eds. 
J. Pochodzalla and Th. Walcher (SIF and Springer- Verlag, Berlin Heidelberg, 2007) p. 243 
|arXiv:nucl-th/0701035| ; W. Weise, R. Hartle, Nucl. Phys. A 804, 173 (2008). 

[20] N.V. Shevchenko, A. Gal, J. Mares, Phys. Rev. Lett. 98, 082301 (2007). 

[21] N.V. Shevchenko, A. Gal, J. Mares, J. Revai, Phys. Rev. C 76, 044004 (2007). 

[22] Y. Ikeda, T. Sato, Phys. Rev. C 76, 035203 (2007). 

[23] E. Friedman, A. Gal, J. Mares, A. Cieply, Phys. Rev. C 60, 024314 (1999). 
[24] M.M. Sharma, M.A. Nagarajan, P. Ring, Phys. Lett. B 312, 377 (1993). 
[25] Y. Sugahara, H. Toki, Nucl. Phys. A 579, 557 (1994). 
[26] Th.A. Rijken, Phys. Rev. C 73, 044007 (2006). 

[27] A. Martinez Torres, K.P. Khemchandani, E. Oset, submitted to Eur. Phys. J. A, 

larXiv:0712.1938 [nucl-th]. 
[28] P. Finelh, N. Kaiser, D. Vretenar, W. Weise, Nucl. Phys. A 770, 1 (2006). 
[29] N.K. Glendenning, J. Schaffner-Bielich, Phys. Rev. Lett. 81, 4564 (1998). 
[30] S.R. Amendolia et al, Phys. Lett. B 178, 435 (1986). 

[31] G. Simon, C. Schmitt, F. Borkowski, V. Walther, Nucl. Phys. A 333, 381 (1980). 



21 



